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O ■ Abstract 

^ . Quantum vacuum energy has been known to have observable consequences since 1948 when 

o : 

I Casimir calculated the force of attraction between parallel uncharged plates, a phenomenon con- 

firmed experimentally with ever increasing precision. Casimir himself suggested that a similar 
attractive self-stress existed for a conducting spherical shell, but Boyer obtained a repulsive stress. 
Other geometries and higher dimensions have been considered over the years. Local effects, and 
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Q I divergences associated with surfaces and edges were studied by several authors. Quite recently, 

O ■ 

Graham et al. have re-examined such calculations, using conventional techniques of perturbative 

(N . 

O ■ quantum field theory to remove divergences, and have suggested that previous self-stress results 



Oh! 



may be suspect. Here we show that the examples considered in their work are misleading; in 
particular, it is well-known that in two space dimensions a circular boundary has a divergence in 
the Casimir energy for massless fields, while for general spatial dimension D not equal to an even 



^ I integer the corresponding Casimir energy arising from massless fields interior and exterior to a 

' hyperspherical shell is finite. It has also long been recognized that the Casimir energy for massive 

fields is divergent for curved boundaries. These conclusions are reinforced by a calculation of the 
relevant leading Feynman diagram in D and in three dimensions. There is therefore no doubt of 
the validity of the conventional finite Casimir calculations. 
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I. INTRODUCTION 



The Casimir effect remains one of the least intuitive consequences of quantum field theory, 
and stands rather outside the usual development of renormalization theory. This is because 
it is inherently nonperturbative, in that macroscopic boundary conditions or backgrounds 
cannot be easily mimicked by perturbative interactions. Its origins go back to the very 
beginnings of quantum mechanics, because it can be thought of as the change in the zero- 
point energy when the background is introduced. 

After examining the van der Waals interaction between two molecules and between a 
molecule and a conducting plate |]T|, Casimir was challenged by Bohr ^] to interpret 
this interaction in terms of zero-point energy [Q, and then to recognize that the zero-point 
fluctuations of the electromagnetic field implied a force between two such plates The 
attractive nature of this force was obviously consistent with the action-at-a-distance inter- 
pretation of it as due to the attraction between fluctuating dipoles making up the material 
of the plates. But intuition flew out the window when Boyer discovered that the energy, and 
hence the self-stress, on a perfectly conducting spherical shell of zero thickness was positive 
or repulsive ||^. Later, it was found that a cylinder was intermediate, giving rise to a small 
but attractive force 0. 

Dimensional dependence was also dramatic. Sen examined a circular boundary in two 
dimensions and found that the energy was infinite^ [^, |T2|. This was later found to be part 
of a pattern: For a hyperspherical shell in D spatial dimensions, the Casimir energy of a 
massless scalar field was finite except when D was a positive even integer, where the energy 
or stress exhibits a simple pole [l^. (For -D < 0, branch points occur at the integers.) An 
intuitive explanation of this, and of the corresponding sign changes, is still lacking. 

Deutsch and Candelas were the first to examine the local effects of fluctuating fields [jl4| , 



for other than the geometry of parallel planes, which was considered by Brown and Maclay 
a decade earlier [jl5|. Typically, surface divergences occur near boundaries, although for 



flat boundaries with conformally-coupled fields, those divergences disappear. The reason 
the global Casimir energy of a (hyper) sphere is finite is that there is a perfect cancellation 
between the interior and exterior divergences. This perfect cancellation is spoiled if the shell 



^ Unfortunately, the author had apparently forgotten this divergence in Ref. [||, wherein an attempt was 
made to extract a finite Casimir energy for a circular boundary. The error was pointed out in Ref. |lOj . 
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has finite thickness, or if the speed of hght is different on the two sides of the boundary [jT6 



Giving the fluctuating field a mass also yields an unremovable divergence [|1^, |1^, [T^, 
except for the case of plane boundaries. 

Recently, Graham et al. ^ have questioned these findings. They have developed 
an approach in which idealized boundary conditions are replaced with interactions with an 
external (nondynamical) field. Potentially divergent terms are subtracted and replaced by 
perturbatively calculable Feynman diagrams. After renormalization of these diagrams, the 
limiting case when the external field becomes a delta function is taken. In this way the 
results for the D = 1 force are reproduced; but the authors find those finite results rather 
unsatisfactory, so they discuss how their limiting procedure gives rise to a different energy, 
corresponding, however, to the conventional force. Then they turn to a circular boundary 
in two space dimensions and find that the Casimir energy is divergent; the implication is 
that this is a general feature, so that all calculations of Casimir self-stress are called into 
question. 

However, as we remarked above, D = 2 is a singular point. What is called for is a cal- 
culation for general D. That is the purpose of this paper. For simplicity, our attention will 
be restricted to scalar fields. We will first, in Sec. re-examine the D = 1 calculation, 
and show that the force is completely finite, while the energy density, or more generally, the 
stress tensor, has a constant divergent part which would be present if the boundaries were not 
present, and is therefore quite without observable consequence. For general D, unphysical 
surface divergences appear in the stress tensor (unphysical because they do not contribute 
to the stress on the sphere), which, for zero mass, vanish if the conformal stress tensor is 
used. Then, in Sec. pT| , we re-examine the self-stress on a sphere in three dimensions, using 
time-splitting to regulate the divergences. The result is, once again, unambiguously finite. 
The critical calculation is given in Sec. |V|, where we review and simplify the diagrammatic 
subtraction method, and explicitly compute the graph in which two external fields are in- 
serted, in D spatial dimensions. As expected, the result is divergent at D = 2,4,6,..., 
but is otherwise finite for D > 3/2. We verify the finiteness for three dimensions by di- 
rectly calculate the oversubtracted graph for D = 3, where we see that the divergences in 
E are independent of the radius of the sphere. Concluding remarks are offered in Sec. |V]. 
Some discussion of how dimensional continuation extracts the part of Feynman graphs which 
contains the dependence on physically relevant parameters is given in the Appendix. 
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II. CASIMIR EFFECT FOR DIRICHLET PLATES 



A. Massless Scalar in 1+1 Dimensions 

We begin by reconsidering the Casimir effect for a massive scalar field which vanishes 
on two parallel plates (Dirichlet boundary conditions). Although these considerations are 



familiar, and are given in some detail in Ref. p4l, we will concentrate on the local effect in 



1+1 dimensions in order to make the divergence structure manifest and make contact with 
the work of Graham et al. [E3 . 



For a massless scalar field 0, the stress tensor is 

T^"" = d^'cpdy - ^g^^dxcpd^cp. (2.1) 

It will be noted that for one spatial dimension, this canonical tensor coincides with the 
conformal one, 

T% = 0. (2.2) 
The scalar field satisfies the free equation 

-d'^(f) = 0, (2.3) 

but is subject to the Dirichlet boundary conditions on the plates at a; = and x = a: 

(j)(x = 0) = (pix = a) = 0. (2.4) 

The corresponding Green's function satisfies 

-d^G{x, t; x', t') = 6{x - x')6{t - t'), (2.5) 

and 

G(0, t; x', t') = G{a, t; x', t') = 0. (2.6) 

Since the Green's function is translationally invariant in time, it is natural to introduce a 
corresponding Fourier transform, 

Gix, x'- t-t')= r (*-*')^7(x, x'; u)- (2.7) 



oo 



the reduced Green's function satisfies the ordinary differential equation 

- (uj^ + ^) gix, x'- u) = 5{x - x'). (2.8) 



We only need the solutions of this equation in two regions: 

, / / \ sinc<jx< sinc<j(x> — a) 

0<x,x < a : g[x,x;u!) = , (2.9a) 

u sm ua 

a<x,x': ^(x,x';tu) = -sinw(x<-a)e^l'^l('^>-''\ (2.9b) 

to 

Here a;> (x<) is the greater (lesser) of x and x'. These are to be compared to the free Green's 
function, when no plates are present: 

g,{x,x';u;) = -^e'\-\^'-^\ (2.10) 

When we recognize that the Green's function is the time-ordered product of the fields, 

(0(x, t)(p{x', t')) = -G{x, t; x', t') (2.11) 

we see that the vacuum expectation value of the stress tensor may be obtained by applying 
a differential operator to the Green's function, and then taking the spacetime points to be 
coincident. For the 00 component, that is, the energy, that differential operator is 

dA + \d'^d'^ = \do&, + \d.d'^. (2-12) 



and so we obtain between the plates 



smciJxsmco'(x — a) + coscija; cosco'(x — a)\ 

2tc 2 usinuja 

duj ioj 

cot uja 

27r 2 

1 /"^ 

^ ~4^y f^CCcothCa. (2.13) 

Here, in the last step we have made the complex frequency rotation, 

u^iC (2.14) 

We notice that this last integral in Eq. ( p.l3| ) does not exist. This is because for large C, 
the hyperbolic cotangent approaches unity. If we subtract off this limiting value we obtain 
a finite result: 



1 1"^ 

-— rfCC(cothCa-l) 



1 1 



o2C,a 



The energy is obtained from this by multiplying by the distance between the plates: 

E = - — , (2.16) 
24a' ^ ' 



which is the well-known Lixscher potential p5[] . 

In the same way we can calculate the vacuum expectation value of the xx component of 
the stress. The relevant differential operator 

d.d'^-\d,d"^\{u' + d.d'^) (2.17) 

is unchanged, so we obtain the same result as for (T°°). The off-diagonal terms (T°^) result 
from the application of the symmetric differential operator 

i(9°9'^ + 9^9'°), (2.18) 

so are necessarily zero. Keeping the divergent term we subtracted off, the result for the 
stress tensor between the plates, < x, a;' < a, is 



TT 



2Aa? 



1 



(^n= ^vac-TT-. I h (2.19) 
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where the divergent terms is 

Mvac = -;^ / dCC- (2.20) 
2vr Jo 

Note that (T^'^) is traceless, 

(T%) = 0, (2.21) 

as required by conformal symmetry. 

If we follow the same operations to find the stress tensor outside the plates from Eq. ( |2.9b|) 
we obtain 

(T00) = (T,J = - ! — - r2cu|cu|coscu(x-a)e'l"l("-") +u;2gj^^/^_^Ng*|..|{.-a)i 

2i J 27C to -' 

= — dcu\u\ = -— c/CC = Wvac. (2.22) 
That is, in the two regions outside the plates, x < or x > a, 

{Tn = M^ac I ^ ° I . (2.23) 



This is exactly the stress tensor that would be found everywhere if the free Green's function 
go in Eq. ( |2.1CI| ) were used. This means that the force on one of the plates is completely 
finite and unambiguous, because it is given by the discontinuity of the xx component of the 
stress tensor across the plate (which follows immediately from the physical meaning of the 
stress tensor in terms of the flux of momentum): 

' TT 



x=a+ 



24a2' 



(2.24) 



Since energies are undefined up to a constant, without any loss of generality we may take 
the stress tensor to be completely finite: 



(T^-^ix)) 



1 



24a2 



1 
1 



< X < a, 
X < or X > a. 



(2.25) 



B. Massless Scalar in 3+1 Dimensions 



In higher dimensions, surface divergences appear. These were discussed in detail in 



Ref. [^, §11.1, but for the sake of completeness we repeat the discussion here. 

In three space dimensions, the use of the canonical stress tensor ( |2.1| ) leads to the following 
expression for the vacuum expectation value of the energy density, 

duj (Pk 



(2.26) 



2n (27r) 

where we have Fourier transformed both in frequency and transverse momentum. If we take 
the plates to be located at ^ = and at 2; = a, we obtain by applying the differential 
operator 



to the Green's function ( 2.9a ) with uj ^ X = yZP^—W^ 

sin Ax< sin A(x> — a) 



< X, x' < a : g{x, x'; A) 



A sin Aa 



The result, 



1 



2%\ sin Aa 

is evaluated by making a Euclidean rotation 



cos Aa — A;^ cos A(22 — a)]. 



(2.27) 



(2.28) 



(2.29) 



a;^<, A 



(2.30) 



and introducing polar coordinates in the (, k plane, 



so 



( = K,cos6, k = KsinO, (2-31) 



(T°°)(z) = I KdK dOn'' " [cos' g cosh /ta 



47r^ Jq Jq sinh na 

+ sin' 9 cosh k{2z — a)] 

1 f°° 1 

/ dn [cosh Ka + 2 cosh fi;(22; — a)l 

127r2 Jo sinh/ta^ ^ 

/ + ^ + ^ 7^ ^ . (2.32) 



67r2 Jq ye^^" - 1 2 e^^" - 1 

Notice that the second term in the last integrand here corresponds to a constant energy 
density, independent of a, so as before it may be discarded as irrelevant. If we integrate the 
third term over 

^ dz [e"^^ + e"'^''-')] = [e"^'^ - l] , (2.33) 

we obtain another (divergent) constant term, so the only part of the vacuum energy corre- 
sponding to an observable force is that coming from the first term: 



00 

^" ^0 



dz LT'''){z) = / dn- = (2.34) 

^ ' 67r2 L e2'^« - 1 1440a3' ^ ' 



which is the well-known Casimir energy/area for a massless scalar field subject to Dirichlet 
boundary conditions, one-half that for an electromagnetic field [|^. 
In general, we have 



(rO°)(^) = u + g{z), (2.35a) 



where 



vr' 



u = T, (2.35b) 

1440a4' ^ ' 

1 1 f°° ^^y^/a ^y{l-z/a) 



giz) = / dyy" ' . (2.35c) 

' 67r2l6a4 ev-l ^ ' 



If we expand the denominator in a geometric series 

1 



ey -1 1 - e-y 

n=l 



5^e-"^ (2.36) 
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FIG. 1: The singular part of the local energy density between parallel plates at z = and z = a. 



we can express g in terms of the generalized or Hurwitz zeta function, 

C(s, a) = 7 — ■ — -, a 7^ a negative integer, 

as follows: 



n=0 



■[C(4,z/a) + C(4,l-;^/a)]. 



(2.37) 



(2.38) 



This function is plotted in Fig. |1|, where it will be observed that it diverges quartically as 
2; — 0, a. (Its z integral over the region between the plates diverges cubically.) As we have 
seen, this badly behaved function does not contribute to the force on the plates. 

Next, we turn to (T^^). According to the stress tensor ( p^.l[ ) and the Green's function 
), that is given by 



duj (fk 



{OA + A 



2i J (27r)3 

1 f duj(fk A 



sin Xz^ sin X{z^ 



2i J (27r)'^ sinAa 



A sin Xa 

[cos Xz cos X{z — a) + sin Xz sin X{z — a) 



(2.39) 



which is independent of z; that is, the normal- normal component of the expectation value of 
the stress tensor between the plates is constant. If once again, the irrelevant a- independent 



9 



part is removed,^ what is left is just three times the constant part of the energy density 

(Em), 

2 

(T,,) = -3 X 



1440a4 

The remaining nonzero components of the stress tensor are 



(2.41) 



2i 

1 f duj <fk 1 



2i J (27r)3 AsinAa 



[uj"^ sin Xz sin A(z — a) 
- A^ cos Xz cos A(2; — a)] 



-u-g{z), 



(2.42) 



where we have again introduced polar coordinates in the frequency-wavenumber plane, and 
again dropped the infinite (a-independent) constant in u. Thus the tensor structure of stress 
tensor is 



(2.43) 



/l o\ 
0-100 
0-10 
\0 3/ 

where u is given by ( p.35b| ) and g by p.38| ). Because u is constant, this vacuum expectation 
value is divergenceless, since g{z) does not contribute to (T^^): 



/l o\ 
0-100 
0-10 
\0 0/ 



(2.44) 



The second term in (|2.43|) diverges at the boundaries, z = 0, a, and has a integral over the 
volume which diverges; yet as we have seen, it is physically irrelevant because its integral 
is independent of a, and it has no normal component. Is there a natural way in which it 
simply does not appear in the local formulation? 

The affirmative answer hinges on the ambiguity in defining the stress tensor.^ It was 



noted in Ref. [24 that this ambiguity was without effect as far as the total stress or the 



The infinite parts of (Tqo) and (Tzz) are related by the same factor of three as the finite parts: 



\inf 



47r2 



, (Too)' 



dK 

127r2 



(2.40) 



For a rather complete discussion of this see Ref. p9], Sees. 3-7, 3-17. 
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total energy was concerned. Now, however, we see the virtue of the conformal stress tensor 

f^"" = d^^dy - -g^"'dx(pd^^ - -(d^'d^ - g^''d^)(p\ (2.45) 

6 



which, because of the equation of motion = 0, has a vanishing trace. 



T% = 0. 



(2.46) 



If we use this stress tensor rather than the canonical one, we merely need supplement the 
above computations by that of the vacuum expectation value of the extra term. Thus to 
obtain (f ^•^■) we add to (IjD 



-{d'^ + dl-d',)G{x,x) 



1 

6i 



duj d?k 



(27r)3 [ AsinAa 
duo d^k 2A 
6z J {2'kY sinAz 



sin \z sin A(z — a) 



cos \{2z — a) 



(2.47) 



which just cancels the surface- divergent term in ( |2.42| ). Again, because G(x, x) only depends 
on z, there is no extra contribution to {Tzz)'- 



-\{dl - g..d'W) = lidl + dl - dl)G{x,x) = 0. 

The extra term for (Too) is just the negative of that in (|2.47| ), 

-^^dlG{x,x) = -g{z), 



(2.48) 



(2.49) 



(2.50) 



which cancels the second term in ( p. 35a] ). Thus, the conformal stress tensor has the following 
vacuum expectation value for the region between the parallel plates: 

/l o\ 
0-100 
0-10 
\0 3/ 

which is traceless, thereby respecting the conformal invariance of the massless theory. This 
is just the result found by Brown and Maclay by general considerations |15[, who argued 
that 

(f'^^) = M[4^''r -^'^^l, (2.51) 
where is the unit vector in the z direction. 
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C. Massive Scalar in D Spatial Dimensions 



It is instructive to repeat the above calculation for a massive scalar field where the plates 
have D — 1 transverse dimensions. We will use the conformal stress tensor, 

T"'' = a'^^a^ - ^gi'''{dx(j)d^(f) + f?^^) - aid'^d" - g''^d^)(l)\ (2.52) 

Here a has to be chosen to be {D — l)/(4i5) in order that the trace vanish (by virtue of the 
field equations) in the massless limit: 

« = ^: T% = -/iV. (2.53) 

The calculation proceeds very similarly to that given above. The only new element is writing 
the momentum integral in polar coordinates: 

dD-^k= \^ fc^'^cifc, (2.54) 

r(i^) ' 



and then introducing polar coordinates as in Eq. (|2.31|) . We encounter the integrals 



relative to which 



rf^(sin^)^- = 2^-^i^L^, (2.55) 



(sin2^) = ^^, {co^H) = ^. (2.56) 



The result for the various nonzero components of the stress tensor are (k^ = + yU^) 

(TOO) ^ - / dpp^-' . [p^cosh/.a + /i^cosh/.(2z-a)], (2.57a) 

iJl (-L'/2) Jq Ksmh/ta 

2-D^-D/2 r-oo 



(T^^) = {Tyy) = ... = -(T°°). (2.57c) 

Surface divergent terms, which do not contribute to the observable force, appear proportional 
to the square of the mass. Of course, the trace of the expectation value of the stress tensor 
is nonzero because of the mass: 

2-D^-D/2 r-oo ^ 

= —ji^ — - — -— / dp p ^ [cosh Ka — cosh k(22; — a)]. (2.58) 

V[D/2) Jq Ksinh/ta 
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As before, the infinite a-independent stress which would be present if the boundary were not 
present (coth/ta 1) may be removed, as it does not contribute to the force on the plates. 
The well known expressions for the force and the energy between parallel plates may be 
easily recovered. We do not see the necessity for the additional terms found by Graham et 
al. [^] in the energy to make the energy finite at zero separation (the fact that the Casimir 
energy diverges at a = reflects the infinite amount of energy released when the plates 
are pushed into coincidence) nor the requirement that the energy should be infinite at zero 
mass, when the observable force is finite there. 



III. SCALAR CASIMIR EFFECT FOR A DIRICHLET SPHERE 



The calculation given in Sec. 11 A was that for a sphere in one spatial dimension. Now 



we consider a massless scalar in three space dimensions, with a spherical boundary on which 
the field vanishes. This corresponds to the TE modes for the electrodynamic situation first 
solved by Boyer 0. The general calculation in D dimensions was given in Ref. [T^; the 
force per unit area is given by the formula 

^(2l + D- 2)T(l + D-2) r°° , , r. / X / X 9 m / x 

Here u = I — 1 -\- D/2. For D = 3 this expression reduces to 

1 °° d 
•^ = -^;^ 5^(2^ + 1) / dxx—ln[li+y2{x)Ki+y2ix)/x]. (3.2) 



In Ref. [|13[ we evaluated this expression by continuing in D from a region where both the 
sum and integrals existed. In that way, a completely finite result was found for all positive 
D not equal to an even integer. 

Here we will adopt a perhaps more physical approach, that of allowing the time- 
coordinates in the underlying Green's function to approach each other, as described in 
Ref. ||2^. That is, we recognize that the x integration above is actually a (dimensionless) 
frequency integral, and therefore we should replace 

1 /"°° 

dxf{x) = - dye'y'fi\y\), (3.3) 

^ J-oo 

where at the end we are to take 6-^0. Immediately, we can replace the x~^ inside the 
logarithm in Eq. ( p.2| ) by x, which makes the integrals converge, because the difference is 

13 



proportional to a delta function in the time separation, a contact term without physical 
significance. 

To proceed, we use the uniform asymptotic expansions for the modified Bessel functions, 
as described in detail in Ref. [El]. This is an expansion in inverse powers of u = 1 + 1/2, low 



terms in which turn out to be remarkably accurate even for modest /. The leading terms in 
this expansion are 

zt 1 1 

In [xIi+i/2{x)Ki+y2{x)] ~ In y + —g{t) + —h{t) + ..., (3.4) 
where x = vz and t = (1 + z'^)'^^'^. Here 

git) = l{t'-6t' + 5t'), (3.5a) 

h(t) = ^(13t^- 284^6 + 1062t^- 1356ti° + 565t^2). (3.5b) 
64 

The leading term in the force/area is therefore 

-j^ POO 

1=0 -^^ 

= F = — ^C(-2) = 0. (3.6) 

1=0 

where in the last step we have used a formal zeta function evaluation."^ Here the rigorous 
way to argue is to recall the presence of the point-splitting factor e*^^'' and to carry out the 
sum on / using 

°° . 11 



2i sin z6/2 ' 

1=0 ' 



so 



u e 



d{z6y2smz6/2 



8 V sm^z6/2 sin z6/ 2 
Then jFg is given by the divergent expression 



' +^1. (3.8) 



TT^a^d'^ J_^ Z'^ 1 + Z'' 



^ Note that the corresponding TE contribution for the electromagnetic Casimir effect would not be zero, 
for there the sum starts from ^ = 1. 
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which we argue is zero because the integrand is odd. 

The next term in the uniform asymptotic expansion ( p.4| ), that involving hkewise gives 
zero pressure, as intimated by the formal zeta function identity, 



oo 



5^z.^ = (2--l)C(-s), (3.10) 

1=0 

which vanishes at s = 0. The same conclusion follows from point splitting, as we can see 
through use of the Euler-Maclaurin sum formula, 

E/(0=/ dlf{l) + -J{Q)-Y,T^,f''-'\^)- (3.11) 

1=0 fc=l ^ 

Here we have 

iz5/2 I I 

We argue again that the first term here gives no contribution to the integral over z because 
it is odd, and then the first two terms in the Euler-Maclaurin formula give 

1 



^1 



1 , d , ^ 1 , d , , 



0. (3.13) 



Derivatives of e^^^^ with respect to I all vanish at z = 0. Again, this cancellation does not 
occur in the electromagnetic case because there the sum starts at / = 1. 

So here the leading term which survives is that of order z/~^ in Eq. (^.4[), namely 



1 °° r-oo 

1=0 -^^ 

where we have now dropped the point splitting factor because this expression is completely 
convergent. The integral over z is 

= ^ (3.15) 

and the sum over / is 3^(2) = 71^/2, so the leading contribution to the stress on the sphere 
is 

357r2 0.00527094 

^2 = 4:7Ta^J^2 = n = ^ . 3.16 

65536a2 ^ ^ 

Numerically this is a terrible approximation. 

What we must do now is return to the full expression and add and subtract the leading 
asymptotic terms. This gives 

^ oo 

S = S2-—^J2(2l + l)Ri, (3.17) 

1=0 
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where 



where the integral 



Ri = Qi + 



Qi 



1 1 

dx \nzt + —g{t) + — 



dx \n}(lxIy{x)Ki,{x)] 



(3.18) 



(3.19) 



was given the asymptotic form in Ref. [13 



Qi 



Z/TT TT 



357r 



5657r 



2 128z/ 32768Z/3 1048577//^ 
120 8 7 677r 138 8 3 577r 



/ > 1. (3.20) 
2147483648Z/7 137438953472z/9' ^ ' 

The first two terms in Eq. (|3.20|) cancel the second and third terms in Eq. ( p.l8| ), of course. 
The third term in Eq. ( |3.20|) corresponds to h{t), so the last three terms displayed in 
Eq. (p.20|) give the asymptotic behavior of the remainder, which we call w{v). Then we 
have, approximately. 



^ n ^ oo 

S^S2 2^1^^^ 2 yi VW{V). 

/=0 i=n+l 



(3.21) 



For n = 1 this gives S ~ 0.00285278/a^, and for larger n this rapidly approaches the value 



first given in Ref. |T3 



S = 0.002817/0^ 



a value much smaller than the famous electromagnetic result ^ ^ ^ 

0.04618 



S 



EM 



(3.22) 



(3.23) 



because of the cancellation of the leading terms noted above. 



IV. DIAGRAMMATIC DIVERGENCE STRUCTURE 

In the two previous sections we have come to rather different conclusions from those of 
Ref. . For the case of parallel plates, studied in Sec. ^ we found: 

• The massless theory is perfectly well defined (no infrared divergences), and surface 
divergences, which in any case have no physical consequences, do not appear if the 
conformal stress tensor is used. 
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• The vacuum expectation value of the stress tensor for the case of a massive scalar does 
have surface divergences, which are proportional to the mass squared, but which do 
not contribute to the force and are therefore physically irrelevant. 

For a massless scalar with a spherical boundary in three dimensions, the formal expressions 
for the force/area and the energy are formally divergent, yet if they are regulated, say 
by point-splitting, the divergences cancel and the energy and self-stress on the sphere are 
completely finite and unambiguous. 

The authors in Ref. [53] came to different conclusions. However, their disagreement with 



us on the D = 1 case seems entirely semantic, and without observable consequence. Their 
substantial argument hinged on their D = 2 calculation. However, it is well known that 
the Casimir effect for a circle is divergent, so it is hard to draw general inferences from an 
examination of that situation. Here, we will re-examine some of the general arguments of 



Ref. for a hypersphere in D space dimensions. 

The general analysis for that case was given in Ref. [0; it is clear that the point-splitting 
method given in the previous section could be applied in that calculation. Instead, we will 
here focus on the issue of the second-order Feynman graph which supposedly is the signal for 
the divergence of the theory in any number of space dimensions. (It is the oversubtracted 
graph which leaves the mode sum more convergent.) We will adopt a somewhat simpler 



formalism than that given in Ref. p3[, based on the "trace- log" formula for the energy. 



E = ^TrlnG, (4.1a) 
where for a "polarization" operator H 

G = Go(l + HG) = Go(l + HGo + UGqUGo + ...). (4.1b) 



The highly sensible approach of Graham et al. ^ is to replace ideal boundary 
conditions by an interaction with an external field a. The Lagrangian for the scalar field is 
thus taken to be 

^ = -lid^4>d^(P + m^^ + a(r)02), (4.2) 

where, anticipating spherical symmetry, we have taken the external field to depend only on 
the spatial radial coordinate. In the end, we may take a to be a delta function, 

a(r) = ^(5(r-a), (4.3) 
a 
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where g is dimensionless and the formal g ^ oo hmit corresponds to the situation of a 
Dirichlet spherical shell. We can now evaluate the one-loop vacuum energy by the replace- 
ment n ^ cr in Eqs. ( [4.1a] ), ( |4.1b| ). It is the second-order graph that is supposed to signal 
nonrenormalizability. 



General D 



We first carry out the calculation in D dimensions. 

E = -—TiaGoaGo 
2T 



I 

2f 



d^+^x d^+^y a{x)Go{x - y)a{y)Go{y - x) 



■ f ,D n n n f f d^'p d^'q e*(p-q)-(x-y) 
= m d xd y a[\x\)a(\y\) — 77r-TT7 77r-TT7^-;^ ^wr-^ (4-4) 

where in the last line we have carried out the integral on t and t', and as a result p^ = = uj. 
Now we introduce polar coordinates, so in terms of the last angle 

d^x = Ao^ix^'^dx sin^-2 (4 5) 

where = 27r"/^/r(n/2) is the surface area of a sphere in n dimensions. Then we encounter 
a Bessel function 

^ rf^sin^-2^e^lP-^l^=°'^^ = (^^— ^-^^ ^ v^r (^^^^^^ JD/2-i(|p-q|x). (4.6) 
Thus the Fourier transform of the field cdxl) is defined by 

a{k) = j d^xe'^ ^'aix) 

= k { — \ j dxx^'^JD/2^i{kx)a{x). (4.7) 



(This agrees with the expression in Ref. p3[ for D = 2.) 
The expression for the energy reduces to 

' ^(|p-q|)'- (4.8) 



, duj fd^qd^p 1 

E = ZTC — 



27r J (27r)2^ {p^ + m^) {q'^ + m'^) 
We carry out the momentum integrations by first using the proper-time representation to 
combine the denominators: 

ds ^_g'g-s(p^+m.2)-s'(g2+m2) 

Jo 

oo pi 

dsS ^^^-sm^-sil~u)p^-suq^^ (4g) 

Jo 
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+ q^ -|- m 



where in the second hne we replace s — s(l — u), s' su. In terms of k = p — q, we 
complete the square in the exponent by writing 

s(l - M)p2 + sMq2 = s[(p - vkf + \i^u{l - m)], (4.10) 

while the corresponding components combine to give —su'^. Now the frequency and p 
integrals are just Gaussian: 

ldu;e^-'=z^, Jd^{p-uk)e-^^^--^^' ={^y^\ (4.11) 

Finally, we introduce polar coordinates for the k integration, with the result 

E = -2-^D^-D+i/2 ^ { 2 ) / dkk^'-'aikf du[m^ + uil-u)kY/^~''/^ (4.12) 
r (y) Jo Jo 

which yields the D = 2 result given in Ref. [p^ . 
If we choose a delta- function potential, 

a{x) = -5{x-a) (4.13) 
a 

we obtain 

E = rdaJl/2-Ai) tdu[m'a' + eu{l-u)r-'y\ (4.14) 

^ \'2) ^ -J^ 

This appears to converge for < < 2 except for the exceptional case m = 0. In that case 
the u integral is simply 

^ \ 2 ) ^ n^-^-ji^l"^ ^ 2 ' 



v{D-i) r(f 



(4.15) 



and the integral over the Bessel functions is 



POO 

/ die-'Ji,2-Ai) = 2^-' 

Jo 



T{2- D)T{D -3/2) 



r(^)^r(i) 
_ r(i-D/2)r(D-3/2) 

which is valid in the region 



(4.16) 



^<D<2. (4.17) 



Thus the energy for a massless scalar is 



E = -2^-^^- ^ ^ ^ ^ ^-^ --, 4.18 

r(f)^ 
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which we take to be the appropriate analytic continuation for all D. This exhibits poles at 
D = 2,4, 6, ... , in congruence with the known divergence structure of the Casimir effect. 
There are also poles occurring at D = 1, —1, —3, . . . , and at D = 3/2, 1/2, —1/2, .... These 
latter two sequences of divergent dimensions correspond to infrared divergences that have 
no counterpart in the Casimir calculations, unlike the ultraviolet, even-integer poles. For 
space dimension between 2 and 4 the Casimir energy is completely finite, in concert with 
this diagnostic. The divergence at D = 2, even putting aside the question of mass, is seen 
not to be generic. 



B. D = 3 

Instead of dimensional continuation, one can work directly in D = 3. Let us regulate the 
theory by inserting a lower limit Sq — > in the proper-time integration, so that for m = 
the energy ( |4.12| ) becomes 



poo rl 

E = / dkk'^d{kf / du\Yi[sQk'^u{l-u)] 

2 Jo Jo 

^ poo d 

= / dk^dikf^ / du\s^k''uii - u)r 

2%'^ Jo da Jo 



(4.19) 



a=0 



If the derivative acts on anything but k we have 



OO POO 

dkk^a{kf = i2nY / dxx^aixf. 
Jo 



(4.20) 



This diverges as a{x) {g/a)6{x — a); but if we regulate the divergence by point-sphtting 



(j{x) lim cr(x — ,^)cr(x + ^), 



(4.21) 



we have 



poo 

/ dA;fc2^(A;)2 = (27r)VW, 
Jo 



(4.22) 



which is seen to be a contact term, independent of a. 
We are left with 



E 



1 d 
2^iT^ da 



dk k^dikYk^"" 



1 /'OO ^ POO 

— / dxxa{x) I dyy a{y)—— / dk k'^°' sin kx sin ky 
2vr Jq Jo Jo 



(4.23) 



«=o 
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Here we have used the Fourier transformation expression ( |4.7] ) , but replaced Bessel functions 
of order 1/2 by the corresponding trigonometric functions. The k integral is now obtained 
from (—1 < a < 0) 

r(a + l)cos(a + l)5 
dx x° cos px = — ^, (4.24) 

so that the energy becomes 

where we have omitted another infinite term that is independent of a. The result is exactly 
the D = 3 value of Eq. (|4.18| ). The justification for omitting (infinite) constant terms in the 
energy is that they are unobservable, not corresponding to a self-stress on the sphere. See 
also the Appendix. 

We further might observe that this energy (f4.25| ) could not be rendered finite were a 
Gaussian profile rather than a delta function employed. This is not surprising. Finiteness is 
only anticipated for an infinitesimal shell, and not for a smooth boundary with continuously 
changing properties. For example, the Casimir energy for a thick dielectric shell apparently 
contain irremovable divergences. 



V. CONCLUSIONS 



The challenge set forth in Ref. and elaborated in Refs. is physically appropri- 
ate and timely given the development of our understanding of the Casimir effect. Certainly 
those authors are justified in objecting to the loose use of the term "renormalization" in con- 
nection with various dubious processes for removing divergences in boundary-value Casimir 
problems. However, it is important to separate the wheat from the chaff. The Casimir force 
between parallel plates, the self-stress (or the force per unit area) on a perfect (Dirichlet 
or Neumann) spherical or cylindrical shell due to massless fields, the energy of fields 
confined to a curved manifold (a hypersphere or torus for example) [^, ^ ^ ^ are ex- 
amples where the Casimir energy is unambiguous and finite, except for exceptional numbers 
of spatial dimensions. 

Of course, these are special cases, and generically Casimir energies are infinite. This is true 
if fields bounded by a spherical shell have mass, if the shell has finite thickness, or if the speed 
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of light inside and outside the shell has different values. The latter case is the interesting 
one of a dielectric ball, first considered in Ref. |T^. The stress or the energy in that case is 



quartically divergent. It was argued, very tentatively in Ref. |T^, and more forcefully later 
|T^, that the divergent terms could be reabsorbed into the definition of physical properties of 
the material medium, the mass density, surface tension, and the like. This "renormalization" 
was in the spirit of the first use of renormalization in physics PB| , P7|] . Obviously, this was 
not a very convincing argument, and was not on a par with perturbative renormalization of 



a quantum field theory.^ However, fairly recently, the discovery by several groups [39, 40, 
ill , |42| , |43|] that the finite part of the Casimir energy for a dilute dielectric sphere was unique, 
and coincided with that obtained by a regulated (dimensionally continued) calculation of the 



van der Waals energy , did provide some evidence that the divergences could be removed 
unambiguously, and had the practical consequence of destroying the hope of explaining 
sonoluminescence on the basis of quantum vacuum energy ^5 . 



Obviously we are still at the early stages of understanding quantum field theory. The 
nature of divergences in vacuum energy calculations is still not understood. However, there 
are a few established peaks that rise above the murky clouds of ignorance, and we should 
not abandon them because the rest is obscure. 



APPENDIX: DIMENSIONAL CONTINUATION OF FEYNMAN DIAGRAMS 

In spite of its impressive successes, one might have concern about the use of dimensional 



continuation to evaluate divergent Feynman diagrams, such as those considered in Sec. |V. 
What is the meaning of the dimensionally continued expression in dimensions where the 
formula gives a finite result, even though the Feynman integral is manifestly divergent? 
Here we give a simple example of what is going on. 

Consider a theory in d spacetime dimensions. The lowest-order self energy diagram 
gives 

S(^) = -UXfi'-'l,, (A.l) 
where by the well known dimensional regularization formula (trivially derived by simple 



The contrary opinion is expressed by Ref. 



22 



proper-time manipulations) 

d'^l 1 m'^-^ „ / d 



r 1 - • (A-2) 



This equality is derived assuming d < 2. The right-hand side of this equation diverges for 
even d > 2. Elementary field theory treatments of scalar field theory have no hesitation in 
accepting that Eq. (|A.2|) is valid in the neighborhood of c? = 4. Therefore, we should ask 
what does it mean for = 3? There, the dimensionally continued formula says 

If we put in a large momentum cutoff A, we can compute the Feynman integral directly: 

/"^,2 1 Am ^, 

^ dl^^^—. = — - —. (A.4) 



' {2tt)^ Jo P + m? 27r2 47r 
Of course, the integral is linearly divergent as A ^ oo, yet the m dependence is correctly 
captured by the continued result ( [A. 3D . The same conclusion is drawn if other regularization 



schemes are employed, such as a proper-time cutoff. On the other hand, dimensional regu- 
larization says J4 is infinite, and indeed, there is no well-defined finite part of that integral: 



(27r)4 7o P + m? 167r2 V m 
More generally, consider 

dH 1 
{2txY (/2 + m 

For a > and d/2 < a we have 



Hd,a)= I , ■ (A.6) 



771 

Iid,a) = ^r{a-d/2). (A.7) 



This assigns the finite value to the integral with a = 2 and d = 5: 



TTl 

7(5,2) = (dimensional continuation), (A. 8a) 



which is divergent when a momentum cutoff is used: 

A m 
12^ ~ 16^ 



A5,2) = ^-— ^, A^cx). (A.8b) 



The general conclusion is that dimensional continuation gives the correct finite part of 
the Feynman graph. (The reader is invited to examine other, more complicated examples.) 
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Since that finite part is, at least in the cases we consider in this paper, the only part that 
contains reference to physical parameters, e.g., the radius of the sphere, we conclude that it 
is effective in isolating the physically observable energy. 
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